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1. Introduction 

Exciton is an elementary excitation of a semiconductor which consists of a pair of two 
correlated fermions, the electron and the hole. Analogous to the Hydrogen atom, it 
is characterized by a binding energy Eb and a Bohr radius a^. Because an exciton is 
composed of two fermions, it is a composite boson. Particularly, in a bulk semiconductor, 
when the excitation of system is dilute, i.e. nexO.% <^ 1, where Hex is the exciton density, 
a bosonic description of system is convenient [Ij. Also Bose- Einstein condensation of 
excitons, which is an essential characteristic of boson systems has been considered 
theoretically [2]. When density of excitons increase, the above condition is violated. 
In this situation the statistics of excitons deviates from Bose statistics. 

In low dimensional semiconductor systems such as quantum well (QW), quantum 
wire and quantum dot (QD), due to the small dimensions and loss of translational 
symmetry, exciton excitation differs from exciton in bulk materials. In semiconductor 
nanostructures the size of the system strongly affects exciton properties. For example, 
in the case of quantum well it is shown [3], if the well width is larger (smaller) than 
the Bohr radius of exciton, the spectrum of quantum well has properties similar to the 
situation in which excitons are boson (fermion). Hence, the size of the system directly 
affects the quantum statistics of excitons in that system. Recently similar results has 
been obtained for QD In Ref. |1] the effects of different statistics of excitons on 
emission spectra of a QD is investigated, and the origin of different statistics of excitons 
is considered. The same results have also been obtained for the quantum well. If the 
size of QD is smaller (larger) than the exciton Bohr radius, excitons behave like fermion 
(boson). Real statistics of excitons in the interaction is considered in [5] and references 
therein. As mentioned before in high density regime, exciton statistics deviates from 
Bose statistics. This is due to the increase of mutual forces between the excitations of 
the system and then the Pauli exclusion principle plays a dominant role [6] . Appearance 
of Bose statistics of exciton-biexciton system and Pauli exclusion effects in superlattice 
has been considered experimentally 

Bosons and fermions are the only two kinds of particles realized in nature. The 
conditions mentioned for excitons (in one regime they are like bosons and in another one 
like fermions) are property of a special kind of statistics called intermediate statistics 
[8]. Bose and Fermi statistics are two limiting cases of this statistics. Properties of 
this statistics have been considered by many authors [9] — [H]. Operator realization 
of intermediate statistics is similar to g-deformed operators [12]. Bosonic g-deformed 
operators [13] are a generalization of the Heisenberg algebra obtained by introducing a 
deformation parameter q. Deviation of this parameter from 1 shows deviation of algebra 
from the Heisenberg algebra. It is shown that it is possible to describe correlated fermion 
pairs with g-deformed bosons [1A\. Therefore it is reasonable to consider an exciton 
system as a g-deformed system. We assume the creation and annihilation operators of 
excitons obey a g-deformed algebra. A g-deformed description of Frenkel exciton has 
been considered recently |15] . 



3 



The algebra generated by g-deformed operators are given by 

Kb% = hbl-Q-'blh = <l^ (1) 
[f^^H] =bl [h,b^] = -bg. 

where h = V^h is the usual particle number operator. Representation of this algebra is 
given in [1^. In the case of excitons, g-parameter can depend on excitation number and 
physical size of system. 

In this paper we consider the interaction of light with a QD embedded in a 
microcavity. By considering excitons in QD as g-deformed bosons (case of g-deformed 
fermion is straightforward) we study the emission spectrum of the system. As is clear, 
the commutator ([1]) explicitly depends on the number of excitons. Hence, this is a 
system in which light interacts with a nonlinear active medium. Therefore, we shall 
obtain the linear and nonlinear response of a g-deformed exciton system. Knowledge 
of interaction of light with a nonlinear medium (g-deformed excitons) and its optical 
response is important for the interpretation of experimental results such as [30]. On 
the other hand, we compared the obtained results with some experimental ones and in 
this manner we investigate the physical origin of g-deformation of excitons. In section 
2 we derive the spectrum of QD when one exciton mode interacts with a single mode 
cavity-field. In section 3 we consider the interaction of two exciton modes with a single 
cavity mode. In section 4 the nonlinear response of QD is derived up to second order 
of approximation. Finally we summarize our conclusions in section 5. 

2. Model Hamiltonian 

We consider a QD embedded in a microcavity which interacts with a single mode cavity- 
field. We assume the excitations in QD have an intermediate statistics [1], and their 
creation and annihilation operators obey g-deformed algebra. We can express the q- 
deformed operators in terms of ordinary boson operators by the following maps 

K = hl4^^y K = J4^^h\ (2) 

\l n{q- q ^) \l n{q - q ^) 

where h and are the ordinary boson operators and h = V^h. Ordinary commutator of 
g-deformed exciton operators is 

[h,:h\] = ^^[q- + q-^-'-'^]^k{n). (3) 

Deviation of this commutator from ordinary boson algebra ([&, fo^] = 1) relates to 
deviation of g-parameter from 1. It is clear, this generalized commutator depends on 
the number of excitations. It seems that by using this algebra we can consider some 
nonlinear phenomena in the system related to the population of excitons. For example, 
biexciton effects can be considered in this manner as an effective approach. So that the 
deformation parameter g can represent some physical parameters such as the ratio of size 
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of system to the Bohr radius of exciton. Interaction of QD with single mode cavity-field 
in rotating wave approximation can be described by the following Hamiltonian 

H = hjja)a + fuxJexPJ^q + hg{ab]j + ci^bq), (4) 

where a and are creation and annihilation operators of cavity field and [Sj, aj] = 6ij. 
We shall consider a phenomenological damping for the system which relates to both 
subsystems: photon and exciton. As is clear from the Hamiltonian (jlj), the exciton 
number is not a constant of motion. Because of the dependence of exciton operator, 
bq, on the exciton number, resulting equations of motion become a nontrivial set of 
coupled equations. On the other hand, since the total number of excitation (exciton 
and photon) is conserved we can diagonalize the Hamiltonian in the subspace of a definite 
excitation. To consider this dynamics we propose an approach based on diagonalization 
of the Hamiltonian by using the polariton transformation This procedure depends 
on some unitary transformations which diagonalize the model Hamiltonian. As is 
usual in this procedure [18], new operators have the same commutation relation as 
the original operators (free operators). Here, there are two distinct sets of operators, 
the cavity mode operators which obey the usual boson commutation relation and exciton 
operators that are g-deformed boson. Therefore, with the presence of these two different 
statistics, mixed operators (polariton operators) do not have specific statistics. They 
can be considered as ordinary boson operators or g-deformed operators. We consider 
both situations and we study the physical results associated with each situation in the 
resonance fluorescence spectrum of QD. 



2.1. Boson polaritons 

In order to solve the dynamical system, we perform the following transformation 

Pk = Ukbq + VkCb. (5) 

Due to the presence of g-deformed operator 6g, we call this transformation a polariton- 
like transformation. As mentioned before, bq depends on the number of excitons 
explicitly and this causes the Hopfield coefficients and will depend on the number 
of excitons. Hence, the transformation ([5]) can be considered as a nonlinear polariton 
transformation. This kind of transformation has been considered recently for the case of 
Bogoliubov transformation [191 [20]. We assume polariton-like operators obey the usual 
boson commutation relations 

[Pk.vi'] = Skk' [Pk,pl] = \uk\'^k{n) + Ivkl"^ = 1, (6) 

where the operator valued function k{h) was introduced by Eq.(l3]). We choose unknown 
coefficients and Vk so that the Hamiltonian becomes diagonal in terms of the 
polariton-like operators 

H = h^rikplpk, (7) 

k 
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where flk is the polariton spectrum and k refers to different polariton branches. By 
taking into account a phenomenological damping for exciton and photon systems 
separately, the unknown parameters satisfy following set of equations 

[ujexk{n)-nk-i'yex]uk + Vkg = 0, Ukgk{fi) + {uj -Vtk-i'-iph)vk = 0.(8) 

In this set of equations, jex and 'jph are the exciton and photon damping constants, 
respectively. From these equations the polariton spectrum can be obtained as 

_ Uexk{h) + U - i{-fex + Iph) 

fc 2 
± -^J[uexk{n) - UJ- il'jex - lph)V + 4:9^k{n)). (9) 

It is apparent that (/-deformed description of excitons causes the splitting between these 
energy eigenvalues be increased in compare to the case of bosonic description of exciton. 
Using the set of equations (IH]) and the polariton spectrum ([9]) we find the coefficients 
for two polariton branches 



^' k{n)[u - 2Qk + 0Jexk{n) - iijex + Iph)] ' 



Vk 




UJexk{n) - i'^ex - ^k 



- 2VLk + ujexk{n) - z(7ex + iph) 
By employing these coefficients all necessary parameters for the polariton Hamiltonian 
are determined. 

Now we can consider the dynamics of polariton operators. The time evolution of 
polariton operators is governed by the polariton Hamiltonian (171) 

Pk = '^IPk^ H] = -iQkPk- (11) 

Let us consider damping effects by taking into account a phenomenological damping 
term and noise operator in the dynamical equations of polariton operators. Hence, the 
time evolution of polariton operator is given by 

Pk = -inkPk-TkPk + Fp^{t), (12) 

where Fp^, (t) is the Langevin noise operator which depends on the reservoir variables and 
Ffc is the damping constant of kth. polariton branch given by Tk = ^:2f+2VL. Correlation 
functions of the noise operators determine physical properties of the system. The 
Langevin noise operator are such that their expectation values (Fx) vanishes, but their 
second order moments do not [21]. They are intimately linked up with the global 
dissipation and in a Markovian environment they take the form 

(FX(t)4(0) = 2r,5(t-0. (13) 

With neglecting the phonon effects by decreasing the temperature, other sources of 
damping like spontaneous recombination of exciton and photon loss are considered as 
Markovian procedures. It follows, on solving Eq. (|T2|) . that 

Pk{t) = Pfc(0)e(-'^^-^'=)* + /* e(-'^'=-^'=)(*-*')Fp, {t')dt'. (14) 
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In this equation we set initial time equal zero. 

The power spectrum of the scattered light for statistical stationary fields is given 
by [22] 

1 r°° - 

S{r,uj) = -Re {E-{r,t)E+{r,t + T))e"^^dT, (15) 
71" Jo 

where are the positive and negative frequency parts of the electric field operator. 
Expressing field operators in terms of creation and annihilation operators we have 

S{r,u) = ^Re H {a\<d)a{r))e'^^ dr. (16) 



TT 

Here, we set t = 0, and A{r) depends on mode function of the cavity-field. 

Now we can express, the field and exciton creation and annihilation operators in 
terms of polariton ones: 

a = vlpi + V2P2, bq = k{h){ulpi + U2P2), (17) 

and at the time t we have 

a{t)=vlpi{t) + vlp2{t). (18) 

Now to calculate the resonance fluorescence spectrum we have to determine the initial 
state of system, we assume at t = 0, the cavity-field is in a coherent state \a), and 
the exciton subsystem in its vacuum state. Under this condition, by using Eq.(IT4l) the 
resonance fluorescence spectrum is obtained as 



S{r,u) = 



TT 



(19) 



In deriving this result we implicitly assume that at t = the noise operator and polariton 
operators are uncorrelated. Fig.(IT]) shows the plot of S{r,u) versus u for different 
values of deformation parameter q. Material parameters are chosen as a; = 1.75 eV, 
Uex = 1.75 eV, g = 200 fieV, jex = 20 fieV, jph = 40 fieV [23], n = 100 and |ap = 9 . 
As is clear when q = 1, spectrum has similar variation as experimental results [23] • This 
figure shows that when q = 1 (nondeformed case) the power spectrum of the fluorescence 
light is a double peak centered at u = Qi and a; = By increasing deviation of q 
from 1, it is apparent from the different plots in this flgure that splitting between two 
peaks increases and the height of one of peaks decreases. This result has been reported in 
resonance fluorescence of excitons when the biexcitonic interaction is taken into account. 
It has been shown ^ that biexcitonic effects are a red shift of the transition frequencies, 
emergence of sidebands due to the switch-on forbidden transitions and asymmetry of the 
emission spectrum. The binding energy of biexciton in QD causes a shift in the spectrum 
of the system. In the present model the splitting of spectrum (Rabi splitting) depends 
on the q'-parameter. Hence, changing this parameter affects the spectrum. Then as a 
one reason of deviation of excitons from ideal Bose system we can consider Coulomb 
interaction between them. On the other hand, g-deformed exciton operators depend on 
the total number of exciton, and biexciton interaction occurs when there are more than 
one exciton. This similarity makes this clue that the g-deformation can be consider as 
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an effective approach to take into account the biexciton effects. As mentioned before, 
the g-parameter can depend on the size of sample. The plotted resonance fluorescence 
spectrum in Fig.([T]) makes clear some differences of optical properties of different size 
QD. For large values of q, compare with 1, spectrum will reduce to one peak. This case 
is a characteristic of the weak coupling regime. 

2.2. Q - deformed polaritons 

In this subsection we assume that the polariton operators are g-deformed operators. 
According to the g-deformed nature of the exciton system we assume the following 
algebra for polariton operators 

[Pk, pI]s = Pkpi - s-^plPk = s^'^ , (20) 

where s denotes the deformation parameter corresponding to the polariton system and 
ftk shows the number operator for A;th polariton branch. Ordinary commutator for these 
operators is 

[Pk,pl] = lukl^m + \vk\'' = ^[s"^- + 5-^"'=+')] = Mihk). (21) 

Using the same approach of the previous subsection we obtain the following set of 
equations for the coefficients of transformation 

[{UexK'f^) - hex - fl'kM{nk)]Uk + VkQ = 0, 

Ukgk{n) + [u- i-fph - Q'kM{nk)]vk = 0. (22) 
From this set of equations we derive the deformed polariton spectrum as 

^, _ ^exfe(n) + U - i{-fex + Iph) 

' " 2M(nfc) 

, ^[uJexK^) -UJ- ijlex - IphW + Vfe(^)) . . 

and the transformation coefficients read as 



^ k{h)[uj~2Vl'^M{nk)+uj,xk{n)-i{-iex + lph)Y ^ ' 



Vk 



M{nk)[uJexk{n) - i-f^x - ^'kM{nk)] 



Lo - 2f]^M(nfc) + uJexk{h) - ii^^x + Iph) ' 

By determining all the variables, polariton Hamiltonian (diagonal Hamiltonian) will 
be determined. By applying the same procedure as before we derive the resonance 
fluorescence spectrum in this case as follows 

. _ A(r)|ap(Kp + hn ^ r. 

Fig. (|2]) shows the plot of S{r,oj) versus oj for different values of polariton deformation 
parameter s. This figure shows that changes of s-parameter (deformation parameter 
of polariton) does not cause any shift in transition frequencies, but causes strengths of 
peaks increase. 



3. Interaction of light with two exciton modes 

We now consider the interaction of one cavity mode with QD when two exciton modes 
are coupled to the field mode. As before, we assume exciton system is expressed by the 
g-deformed operators. The total Hamiltonian of the system under consideration can be 
written as follows 

H = hwa^a + hJ2 ^eJlS + ^9Y1 (^^L + «^^««)- (26) 

j=l,2 i=l,2 

We assume both excitons have the same coupling constant with the cavity mode. We 
solve this system as before by diagonalizing the Hamiltonian. For this purpose we 
perform the following transformation 

Pk = Ukbg, + Xkbq^ + VkCL. (27) 

We consider the situation in which the polariton operators obey the nondeformed Bose 
statistics 

[Pk,p\] = \uk\'^k{ni) + \xk\'^k{h2) + \vi\^ = 1, (28) 

where hi {i = 1,2) represents the number operator for each excitonic mode. As is 
clear in this case there are three polariton branches. Assuming the transformation fl27j) 
diagonalizes the Hamiltonian fl26|) . this polariton Hamiltonian takes the following form 

H = hY,^kpipk, (29) 

k 

where summation is over all polariton branches. The following equation determines the 
polariton spectrum 

(c - Qk)[{d - Qk){uj - i-fph - ^k) - g'^k{hi)] - g^k{h2){d - Qk) = 0, (30) 

where c = Uexik{hi) — i'jexi and d = Uex2k{h2) — i'jex2- By deriving the polariton 
spectrum the transformation parameters are obtained as 

g[{d-nk){uj - i'jph - ^k) - g^k{h2)] , . 

Uk = , [oi) 

g^k{hi) 

Xk = 



Vk 



A 

c - nk)[{d - nk){u - ijph - ^k) - g^k{n2)\ 



A 

where the parameter A is given by 

A = {[g^k{hi) + (c - fifc)2][(rf - nk){io - hph - ^k) - g^k{h2)]^ 

+ g^k\h,)k{h2)f^. (32) 

In this manner, all the parameters which appear in the polariton Hamiltonian are 
determined. By repeating the approach of previous section the resonance fiuorescence 
spectrum of system with different initial conditions can be determined. If we assume 
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at t = the cavity mode is in the coherent state |a) and QD in vacuum state |0), the 
resonance fluorescence spectrum is given by 

To show complex structure (multi-peak structure) of this spectrum Fig. (|3]) presents 
the spectra on a logarithmic scale. For the sake of clarity, we have powered some peaks 
compare to other ones in this figure. In the case of q = 1 (nondeformed exciton) the 
spectrum has two peaks. Increasing the g-parameter causes that splitting between peaks 
be increased and spectrum becomes multi-peaks. Multi-peaks structure in emission of 
exciton such as Mollow triplet was predicted when excitons obey statistics different 
from Bose statistics [31 H] . When, g-parameter is changed, the energy and intensities of 
emission change. Effects of exciton number on absorption spectrum of QD is considered 
. Due to the relation of absorption spectrum and resonance fluorescence, similar result 
is obtain in I2H . 



4. Nonlinear response of excitons in g-deformed regime 

In previous sections we considered some physical results of g-deformed description 
of excitons. The g-deformed description can be served as a nonlinear description of 
excitons. It is well-known that different kinds of nonlinearity in an exciton system 
lead to different orders of nonlinear response of the system [251 126]. Therefore, we try 
to obtain optical response of a driven quantum dot, which its optical excitations are 
considered as g-deformed systems. For this purpose we will calculate the coefficient 
absorption of a QD in this regime. In this section we neglect all damping effects and we 
consider the Hamiltonian of the system as follows 

H = hjja)a + f\hJexb\i>q + hg^ab}^ + a)bq). (34) 

In the electron picture, the induced dipole moment by transition of an electron is 
described by /i = ajfi-c + aja„ [27|. The operator al (a^,) is the creation (annihilation) 
operator for an electron in the valance band (level in the case of QD), and a\ {a^ is 
the creation (annihilation) operator for an electron in the conduction band. Hence, 
creation of an exciton is denoted by ajcit, = 6j. Therefore we can write the dipole 
operator of QD as // = b\ + bq. The macroscopic polarization is expectation value 
of polarization operator. The optical response function represents the reaction of the 
system to an external classic field E{t) coupled to the variables of system [28], i.g., the 
dipole operator. Hence, we consider an external field as a pump source and we treat 
the reaction of QD to it. Then the total Hamiltonian of system is then given by 

H = hhja)a + hcoexblbq + hg{db}j + a)bq) — [d^c ■ E{t)bq + dcv ■ E{t)P^, (35) 

where d^^ denotes the dipole matrix element. The Hamiltonian in the interaction picture 
has the form 



Hint = % 



(36) 
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The observable of interest for the optical response is the time-dependent dipole density 
li{t) = {bg{t)) + h.c. = Trex{bgPexit)) + h.c, where Tr^x means trace over the exciton 
system and pex(t) = Trfp(t), which p(t) is the total time dependent density matrix of 
the system and pex{t) is the time dependent density matrix of exciton system. The total 
time dependent density matrix is given by 

pit) = Uit,to)pito)U-\t,to), (37) 

where U{t, to) = Texp[— I Hint{t')dt'] is the time ordered evolution operator and p{to) 
is the total density matrix of system at initial time. We assume that the quantum field 
and exciton system are both in vacuum state. Therefore, the time dependent density 
matrix of excitons is given by 

PeM = Y.^n\U{tM)mMe.){eMfmU~\t,to)\n), (38) 



where summation is carried on field state and the matrix elements of the time evolution 
operator are in the basis of field states. By using the Feynman disentanglement theorem 
[29] the matrix elements of the time evolution operator U{t,to) can be evaluated. We 
can write Hamiltonian in ( l36l) as Hint = Hi{t) + H2{t), where 



Hi{t) = hg 

mt) = - 



(39) 



As is clear H2{t) depends only on exciton operators. The time evolution operator can 
be written as 



U{t,to)=Texp[-- I {H,{t') + H2{t))dt'] 



to 



h 



fexp[-- / H2{t')dt']exp[-- / Hi{s)ds]. 



to 



h 



(40) 



to 



In this equation we use Feynman notation [29]. These two exponential terms are not 
disentangle from each other. They are correlated and in doing integration, we have to 
take into account ordering of operators. In calculation of matrix element of this operator 
in the basis of field states, second exponential can be considered as a ordinary c-number 
function of t', because it is independent of field operators: 

{i\U{t,to)\j) = (^|^exp[-^ f H,{t')dt']\j)exp[-l f H2{s)ds]. (41) 



h 



to 



h 



to 



On the other hand we consider all the exciton operators in Hi{t) as ordinary c-number 
functions, and we can write 



{t\U{t,to)\3) = {t\Texp[-- I idle'-'' Bit') + dee~'^''B*it'))dt']\j) 



to 



X exp[-- / H2is)ds] 



(42) 



to 



11 



where B{t) is a ordinary function corresponding to exciton operators. As is clear this 
matrix element is a function of exciton operators. The influence of exciton system is 
completely contained in this operator functional and factored term in f H2l) . By using 
Feynman theorem, above matrix element can written as 



{i\exp[-ig [ ale'''''B{t')dt']exp[-ig [ a[,e-''''' B*{t')dt']\j), 

J to J to 

where in this equation a[, = V^^{t)dtV{t), and 



V{t) = exp[-igd^ / B{t')e 



to 



In this manner the density matrix of exciton system takes the following form 



(43) 



(44) 



(45) 



where 



Si{h,bq,bl) 
S2in,bg,bl) 

L{hex) 



-gbq 



Ju}exk{hex)(t-to) 



^-iuj^xk{hex+l)[t~to) 



bib 



+ 



UJexK^ex + 1) 

Q-ioJex[k{hex+l)-k{hsx-l)]it~to) 

'^"^^[W - UexK^^ex - l)][t^ - UexK'f^ex + 1)] 
p-iuJex[k{h^x+'i)-k{h^x))\(t-to) 

[U - Uexkiflex + 2)] [u - iOexkiflex)] ' 



(46) 
(47) 



and 



/(6„ b\) = f exp [ if dt\l, ■ E{t%e 

^ J to 



+ d,, ■ E{t')b\e 



l\t] Mexk(hex)t' 



-iuJexk{hex + 'i-)t' 



(4J 



By expanding the function f(bq, 6j) up to second order in E{t) and using (H5i) we obtain 
the time-dependent dipole density as follows 

n 



j./\Au}ext' Jj-I 



to 



■^\/ fg{n)ho{n)\e 



Jto 

+ i^,\JUn + l)h,{n + 2)!e-^^(°y/,(n + 2)!/,(n + 2) x 

■ E{t')d^^ ■ E{r)d,^ ■ ^(s)e*'^-[^+"-'=("+^)*']dt'rfrds 

V- 

+ ^yZW^o(n)!e-^^(°)^/,(n)!/,(n) x 
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J to 'J tr\ ^ to 

^f,{n + l)ho{n + l)!e-^^«^/,(n + l)!/,(n+l) f f x 
■ E{r)d^^ ■ E{s)d,^ ■ E{t')e-''^^^^''^"'+^'^^''-''>-''Ut'drds 



'to k) •'to 
i 

2h? 



to ''to J to 



- ^s\fUn)Mn + l)!e-^^«^/,(n + l)!/,(n + l) 
4„ ■ E{r)d,, ■ E{s)dc, ■ E{ty-^^^^''+^^^''-'^+''^dt'drds\, (49) 



X 

to J to J to 



where hi(n) = e( and fq{n) = y^^^r^- These equation shows that in 

this conditions second order response function is equal zero. Now we can calculate 
linear and nonlinear electric susceptibility of this exciton system from this equation. 
Generalized linear and nonlinear absorption spectra of this system is shown in figures 
Q-dn]) for different values of g-parameter. In these plots, Is-exciton is considered. In 
these figures we choose h = e = 1, g = 200 fiev and Uex = 1574mef . Fig. ([4]) shows plots 
of linear absorption spectra and Fig. ([6]) shows plots of nonlinear spectra. On the other 
hand, 3-dimension plot of linear absorption coefficients is given in figure ([5]). It is clear 
that changes of ^-parameter strongly affects absorption spectra of the system. These 
figures show in the presence of g-values absorption of probe beam shows a complex 
structure: a multiple-like absorption pattern appears with one strong peak and some 
side bands. Presence of these side bands is a signature of the optical generation of an 
nonlinear exciton (an exciton which expresses with (^-deformed operator). Negative part 
of the absorption spectrum demonstrates gain of the probe beam. Due to the resonance 
interaction of pump with exciton transition, the gain effect comes from the coherent 
energy exchange between the pump and probe beams through the QD nonlinearity. 
The obtained absorption spectra are very similar to experimental results [30]. In Ref. 
[30] absorption spectra of a driven charged QD is derived experimentally. Charged 
QD is a nonlinear medium and is similar to our model. Then It can be consider as a 
experimental test of our model. 



5. Conclusion 



Q-deformed description of excitons in a QD and its physical consequences was 
considered. We showed that increasing the g-parameter will lead to increase of splitting 
between peaks in the spectrum and asymmetry of spectrum. Similar effects were observe 
when biexciton effects taken into account. In experiments of QD it is shown p3] the 
same results are obtained in different temperatures. Then we can associate this physical 
parameter as source of g- deformation. The temperature dependence of emission energy 
of system can be attributed to the change in the refractive index of its active medium 
with temperature. We have derived the optical response of QD with g-deformed exciton. 
As mentioned before g-deformed description of excitons will lead to dependence of optical 
response on q parameter. Hence, due to the wide range of q parameter and its effects 
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on optical response we can consider some parameters like temperature and interaction 
between the excitons which affects the optical response of QD as sources of g-deformation 
of excitons. As mentioned, the relation of quantum statistics of excitons in the QD 
and the size of QD has been considered. Then we can consider the ratio of exciton 
Bohr radius to dimension of system and exciton population as two main sources of q- 
deformation. Q-deformed operator depends on total number of associated particles of 
system. Therefore we can interpret g-deformed operator as an operator which consists 
of effects of other excitations of system implicitly. Then it is reasonable to consider 
this description as an effective description which takes into account some nonlinearity 
in exciton system. As we saw, in the case of interaction of light with two excitons, 
when q = 1 this system showed a two peaks spectrum. While by increasing deviation of 
exciton from Bose statistics, spectrum becomes multi peak. Due to the nonlinear nature 
of g-deformed exciton we showed that different orders of nonlinear response function of 
this system can be calculated. From coincidence of obtained results and experimental 
results, we can conclude that g-deformed description of excitons can be a considerable 
model for excitons. With comparing the obtained results in this paper with experimental 
ones we can investigate the origin of this description of excitons. As pointed out the 
ratio of system dimension to the Bohr radius of exciton is one of the sources of deviation 
of excitons from usual boson. The obtained results are very similar to the effects of the 
exciton-exciton interaction [3], [31] which is relates to exciton population and biexciton 
binding energy. On the other hand, it is shown that [1] exciton density is another 
source of their deviation from ordinary bosons. To sum up we attribute the origin of 
g-deformation of the excitons to their density, their mutual interactions, confinement 
size and other parameters which cause fluctuation of optical response of the system. 
Q-deformed description of an active medium causes that the optical properties of system 
depend on the g-parameter. Then, it is seem that parameters which can affect optical 
properties of the active medium (like refractive index) their effects can be considered by 
this formulation. The g-parameter can be considered as a variation parameter which its 
values can be obtained from comparison of theoretical and experimental results. 
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Figure 1. Plots of 8(10) versus ui. Parameter are choose as = 1.75 eV, 
LOex = 1-75 eV, g = 200 ^leV, jex = 20 ^J.eV, -fph = 40 ^eV, n = l and jap = 9. Solid 
plot corresponds to q = 1, nondeformed case. Dotted one corresponds to q = 1.01, and 
for dash line q is equal 1.015. 
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Figure 2. Plots of S{uj) versus oj. Parameter are choose as w = 1.75 eV, 
LOex = 1.75 eV, g = 200 /xe]/, jex = 20 iieV, jph = 40 fieV, n = 1 and jap = 9. 
In all figures we have q = 1. Solid line corresponds to case s = 1. In dotted one we 
have s = 1.007 and for dash line s = 1.01. 
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Figure 3. Plots of S{uj) versus uj. Parameter are choose as w = 1.75 eV, 

ujexi = 1.75 eV, u>ex2 = 1-77 eV, g = 200 ^eV, jexi = lex-^ = 200 [leV , -fpu = 45 ^leV , 
77-1 = 1, 77,2 = 1 and |a|^ = 9. Dotted line corresponds to nondeformed case 91,92 — 1- 
For solid line qi, 92 = 1-04. In the case of dashed line gi, (72 — 1.08. 
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Figure 4. Plots of spectrum absorption versus w. We consider Is-exciton and 
Parameter are choose as/i = e = l,(7 = 200 /lew and = 1574 mev. Solid 
plot corresponds to nondeformed case q= \. For dotted one q = 1.01 and in dash one 
g = 0.99. 




Figure 5. 3D-Plots of spectrum absorption versus uj and deformation parameter q. 
Physical parameter are the same as Fig.Q. 




Figure 6. Plots of nonlinear spectrum absorption versus uj. We consider Is-exciton 
and Parameters are choose a&h = e — \,g — 200 /ieu and uj^x = 1574 mev Solid plot 
corresponds to nondeformed case 9=1. In dotted plot q — 1.01. In dash plot q = 0.99. 



